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ABSTRACT 

Context. Protostars of intermediate mass provide a bridge between theories of low- and high-mass star formation. Molecular outflows 
emerging from such sources can be used to determine the influence of fragmentation and multiplicity on protostellar evolution through 
the apparent correlation of outflow forces of intermediate mass protostars with the total luminosity instead of the individual luminosity. 
Aims. The aim of this paper is to derive outflow forces from outflows of six intermediate mass protostellar regions and validate 
the apparent correlation between total luminosity and outflow force seen in earlier work, as well as remove uncertainties caused by 
different methodology. 

Methods. By comparing CO 6-5 observations obtained with APEX with non-LTE radiative transfer model predictions, optical depths, 
temperatures, densities of the gas of the molecular outflows are derived. Outflow forces, dynamical timescales and kinetic luminosities 
are subsequently calculated. 

Results. Outflow parameters, including the forces, were derived for all sources. Temperatures in excess of 50 K were found for all 
flows, in line with recent low-mass results. However, comparison with other studies could not corroborate conclusions from earlier 
work on intermediate mass protostars which hypothesized that fragmentation enhances outflow forces in clustered intermediate mass 
star formation. Any enhancement in comparison with the classical relation between outflow force and luminosity can be attributed 
the use of a higher excitation line and improvement in methods; They are in line with results from low-mass protostars using similar 
techniques. 

Conclusions. The role of fragmentation on outflows is an important ingredient to understand clustered star formation and the link 
between low and high-mass star formation. However, detailed information on spatial scales of a few 100 AU, covering all individual 
members is needed to make the necessary progress. 
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1. Introduction 


In the current view of star formation, low-mass (Mstar <3 M©) 
and high-mass (Mstar >8 M©) star formation are described by 
different theories. Low-mass theories build on the assumption 
of singular collapse ( [Shu et ^|1999| ), while various high- 


mass theories focus on clustering and energetic environments 
( jKrumholz et al. [2009 [ [Tan et al.|2014| ). The poorly studied pro- 
tostars/protostellar clusters of intermediate mass (Mstar between 
3 and 8 M©, Lboi between 30 and 5,000 L©) provide an ideal 
testbed for a unified theory of star formation. Such a theory must 
be able to correctly address the observed level of fragmentation 
within intermediate mass protostellar regions, its influence on 
the observed properties and simultaneously reproduce predicted 
populations set by the initial mass function ( |Bate||2009t jOffnerj 
jet al.|20lot [Hansen al.|201^ . 

Multiplicity and fragmentation of intermediate mass proto- 
sta rs is a relatively unexplored area (for the most recent reviews, 
see [Goodwin et al.|2007|[Beltran|2015[ ). Of those studied in de¬ 
tail, very few deeply embedded intermediate mass protostars are 
truly isolated single protostar s; most are tightly packed clusters 
of low-mass protostars (e.g., [van Kempen et al.][2012| ) with at 
times an actual intermediate mass protostar found near the cen¬ 
ter (e.g., [Euente et al.[[2001[ ). Isolated intermediate mass proto¬ 


stars are known, but likely very rare (LI641 S3 MMSl is the 
best candidate, see [van Kempen et al.|2012 l. 

Interferometric observations at submillimeter wavelengths 
are needed to characterize protostellar content of intermediate 
mass protostars, down to the very young, heavily embedded pro¬ 
tostars. Observations with the required spatial resolution have 


been obtained for just a handful of cases (e.g.. 

Euente et al.[2001 

2005 

Teixeira et al.|2007l Euente et al.|2012[ 

van Kempen et al. 

2012 

Carrasco-Gonzalez et al.[2012|). Most studies lack sensitiv- 


ity and/or spatial resolution to separate emission from individual 
protostars. 


A powerful method to directly probe star formation with¬ 
out time-intensive interferometric observations is through the 
study of bipolar jets and molecular outflows. The bipolar jet 
not only allows part of the angular momentum to disperse and 
gravitational collapse to continue ( [Arce et aT||2007| ), it also en¬ 
trains significant amounts of the surrounding gas, thereby cre¬ 
ating the molecular outflow, which is thought to be the main 
source of mechanical feedback onto the parental cloud ( [Naka-[ 
[mura & Li|[2007[ ). Outflows also facilitate radiative feedback at 
larger radii, as radiation is able to escape more readily through 
the lower density outflow cavities than through the protostellar 
envelopes. These feedback effects dramatically affect fragmen- 
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tation and mass accretion rates (|Offner et al.|2010[|Hansen et al.| 

|20T2l >. 


Outflows strengths are quantifled by the outflow force, Fco- 
Deriving outflow forces is difficult, requiring sensitive observa¬ 
tions of the line wings of molecular tracers across multiple tran¬ 
sitions. In practice, CO is used almost exclusively. Recent bench¬ 
marking limits variations between methods to less than an order 
of magnitude ( [van der Marel et al.|2013| ). 


Outflows emerging from intermediate mass protostars have 
been sparsely studied. The most complete study, [Beltran et al.| 
( |2008| referred to as BOS from here on) studied outflow forces of 
a sample of Intermediate mass outflows using a range of meth¬ 
ods and data in comparison with detailed observations of one 
of them, and proposed a (cor)relation between the Fqo and to¬ 
tal Lboi (see their section 6). This is a change from the relation 
between Fqo and individual luminosity, first identified by [Bon-j 
temps et S^ ( |1996| ). Higher mass accretion rates on the driving 


source, set or influenced by the level of fragmentation, was put 
forward as the origin of this effect. However, observational con¬ 
straints and large uncertainties in the sample and the range of 
methods used to derive outflow forces limited validation of this 
hyopthesis. 


Since the publication of BOS, the Atacama Pathfinder Exper¬ 
iment (APEX) and Herschel Space Observatory have enabled 
regular observations of spectrally resolved n iid- and high-/ CO 
emission lineof molecu lar outflows (See van Kempen et al. 


2009a|b[|Eich et al.|2010t[v^ Kempen et al.|2010t|Yildiz et al. 


2010||2012||Kristensen et al.|2013[|San Jose-Garcia et al.|2Q13j ). 
The reliability of outflow force calculations has improved sig¬ 
nificantly by making use of both the larger range of observable 
CO transitions, the increase in excitation temperatures of these 
transitions, and the ability to map outflows completely within 
a reasonable time. One of the main results is the constraint 
on entrained gas temperatures of individual flows to 50 K or 
van Kempen et al.|2009b|a| ), warmer than the classic ally 


higher ___ 

adopted temperatures of ^^30 K (e.g.,|Bachiller et al.|2Q01 and 


many others), but in line with shock model predictions (jHatchellj 
jet al.|1999D . 


In this paper, we present new observations using the 
CHAMP^ instrument mounted on APE^of CO and ^^CO J=6- 
5 emission of outflows emerging from six intermediate mass pro¬ 
tostars, four of which are known to form clusters of low-mass 
sources. [Cl] 2-1 emission is discussed as a complement. The 
goal of this paper is to validate the relation between outflow force 
and total luminosity proposed by BOS by making use of the ad¬ 
vances of mid-/ CO observations and the methodology applied 
to low-mass protostars developed by this group {van Kempenj 
jet al.|2009b|a]|Yildiz et al.|2QT^ . 

Section 2 describes the observations, source sample and data 
reduction strategy. Results are given in Section 3. Analysis is 
presented in section 4, while we discuss the importance of the 
derived physical parameters in Section 5. Conclusions and future 
work are listed in Section 6. 


^ Throughout this paper low-/ CO transitions are defined as having 
/up <4, mid-/ CO transition having 4 < /up < 9 and high-/ CO transi¬ 
tions having /up >9 

^ This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between the 
Max-Planck-Institut fur Radioastronomie, the European Southern Ob¬ 
servatory, and the Onsala Space Observatory. 
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2. Observations 


The dual-frequency CHAMP^ array receiver ( Glisten et al. 2Q0S I 
mounted on APEX, was used to map the CO J=6-5 and ^^CO 
/=6-5 transitions in six intermediate mass protostars. As a com¬ 
plement, observations of the [Cl] ^P 2 -^Pi were obtained. Obser¬ 
vations were carried out between November 2009 and July 2012 
using the On-the-fly (OTE) mode. East Eourier Transform spec¬ 
trometer back-ends were attached to each of 7 pixels for each 
frequency band, providing a spectral resolution better than 0.1 
km s“^ Typical system temperatures ranged between 1300 and 
2000 K for the 690 GHz receivers and 3500 to 5000 K for the 
SIO GHz receivers. Maps of at least 2'5 by 2'5 (CO /=6-5) and 
V by 1' (^^CO and [Cl]) were obtained. Sensitivities across the 
maps and sample varied by factors of 2-4 because of the dif¬ 
ferent atmospheric conditions in combination with elevation of 
the sources. Noise levels increase by a factor of 2 at the map 
edges (the outer 15 "). The average beam efficiency was derived 
to be 0.4S for the 690 GHz array and 0.42 for the SOO GHz ar¬ 
ray. However, there are variations from month to montlj^ Beam 
efficiencies for individual scans were taken as close in time as 
possible to the observation date. 

Erom the system temperature and beam efficiency measure¬ 
ments the total flux uncertainty is assumed to be 20% for the 690 
GHz band and 30% for the SOO GHz band. 


2.1. Source Selection 


The targeted sample consisted of all intermediate mass proto¬ 
stars observable from APEX fro m the WISH key program on 
the Hersch el Space Telescope p] ( van Dishoeck et al. 2011): 
NGC2071 (jCarrasco-Gonzalez et al.||2012l [van Kempen et al. 
|MT^ , L1641 S3 MMSl ([van Kempen et al.|2012| ). Vela IRS 17 
and Vela IRS 19 ( jGiannini et al.||2005| ). The sample was com¬ 
pleted by Serpens SMM 1 (the most massive low-mass proto- 
star included in WISH, and known to be a single pr otostar. See 
[van Kempen et al.]|2009dt [Kristensen et al.||20'l^ and IRAS 
20050-r2720 (BOS). Table lists all relevant properties (total 
luminosity, distance, total mass, Vlsr and estimated number of 
members). If conflicting values were reported in ex isting litera¬ 
ture, values presented in van Dishoeck et al. ( 201l| p are given. 


2.2. Data Reduction 

During the observations, the raw data-streams were immediately 
calibrated using the APEX on-line calibrator, assuming an image 
sideband suppression of 10 dB. and [Cl] observations of 
L1641 S3 MMSl had to be reprocessed with the APEX off-line 
calibration software owing to inaccuracies in the on-line cali¬ 
bration. Afterwards, a full reduction was done using standard 
routines in the CLASS and GREG packages of GILDAS|^ The 
final data product was transformed into large EITS cubes in main 
beam temperature scale with a spectral resolution of 0.1 km s“^ 
In the [Cl] spectrum of Vela IRS 19, an absorption feature at 
~0 km s“^ is present. This is caused by large-scale cloud emis- 

^ see the MPIfR website for more information and distri¬ 
bution of beam efficiency measurements: http://www3.mpifr- 
bonn.mpg.de/div/submmtech/heterodyne/champplus/ 

Water in Star-forming regions with Herschel, see 
http://www.strw.leidenuniv.nl/WISH for more information 
" or BOS in the case on IRAS 20050+2720 

^ GILDAS is a set of (sub-)millimeter radioastronomical applica¬ 
tions (either single-dish or interferometer) developed at IRAM, see 
http: //w w w.iram.fr/IRAMFR/GILDAS 
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Table 1. Sample of studied IM protostars. The Right Ascension and Declination are the previously derived and/or estimated centers of gravity. 


Source 

R.A. 

(hms [J2000]) 

Deck 

(dms [J2000]) 

kbol 

(Lo) 

Dist. 

(pc) 

Mass 

(Mo) 

f^LSR 
(km s“^) 

# 

members 

Ref.^ 

NGC 2071 

05:47:04.4 

-h00:21:49.0 

520/40^ 

422 

30 

9.6 

5 

1,2 

L1641 S3 MMS 1 

05:39:55.9 

-07:30:28.0 

70/250^ 

465 

20.9 

5.3 

1 

1,2,3,4 

Vela IRS 17 

08:46:34.7 

-43:54:30.5 

715 

700 

6.4 

3.9 

3-f3 

5,6,7,8 

Vela IRS 19 

08:48:48.5 

-43:32:29.0 

776 

700 

3.5 

12.2 

3-f3 

5,6,8 

Serpens SMMl 

18:29:49.8 

-f01:15:20.5 

30 

250 

16.1 

8.5 

1 

9,10,11,12,13 

IRAS 20050-^2720 

20 07 05.8 

+21 29:00.0 

280 

700 

17.1 

6.4 

4+ 

14,15,16,17,18,19 


^ References beyond van Dishoeck et al.|(|2011 

- l:|Wilson et al.H2005 

1 , 2:|van Kempen et al.H2012|, 3: 

Stanke et al.H2000|, 4:|Morgan et al.[ 

199ll as FIRSSE 

101, 5|Liseau et al.|(|1992), 6:|Slawson & Reed 

dl988 

, 7|Giannini et al. 

00 

o 

O 

(N 

Wouterloot & Brandy] 

1999), 9:|Kristensen et al. 

(|2010a), 10: 

Goicoechea et al.| 

\2012), ll:|Hogerheijde et al.|(|1999), 12:|White et al.| 

1995), 13:|Kristensen etal.|i 

2012), 14:|FroebrichH2005), 15:|Bachiller et al.|i 

1995), 16:|Zhang et al.|(|2005). 


iz..|vviiiLC CL ai.|i i j.|rviiaLciiacii cl ai.|i iucuiiciimz.' 

17: |Chini et al.H2001) , 18 |Wang et al.l(p091. 19: B08 _ 

^ Second luminosity numbers from |vaii Kenipen et al.lj2012} , equalling the sum of individual luminosities. 
^ Estimated from the infrared results jOiannini et al.|2005t . 


Table 2. Noise levels, integrated and peak intensities at the central po¬ 
sition. Temperatures are in given in main beam temperature units. 


Source 

CO 6-5 



fTMBdV 

rpeak 

RMS^ 


[Kkms-'] 

[K] 

[K] 

NGC 2071 

819.7 

64.0 

0.27 

L1641 S3 MMS 1 

75.0 

8.8 

0.30 

Vela IRS 17 

241.7 

25.6 

0.9 

Vela IRS 19 

91.2 

7.1 

0.78 

IRAS 20050-h2720 

183.4 

16.2 

0.60 

Serpens SMMl 

151.3 

19.8 

0.20 


‘^CO 6-5 



fTMBdV 

rpeak 

RMS* 


[Kkms-'] 

[K] 

[K] 

NGC 2071 

141.1 

22.8 

0.40 

L1641 S3 MMS 1 

10.2 

5.0 

0.60 

Vela IRS 17 

63.1 

14.8 

0.12 

Vela IRS 19 

10.8 

3.6 

0.12 

IRAS 20050+2120 

10.3 

3.4 

0.44 

Serpens SMMl 

25.5 

7.1 

0.21 


[CI] 

^P2-^Pl 



fTuBdV 

rpeak 

RMS* 


[Kkms-'] 

[K] 

[K] 

NGC 2071 

41.8 

4.0 

0.80 

L1641 S3 MMS 1 

- 

- 

1.90 

Vela IRS 17 

41.2 

7.5 

0.33 

Vela IRS 19 

15.2 

3.3 

0.39 

IRAS 20050-h2720 

12.0 

4.3 

1.30 

Serpens SMMl 

10.2 

2.9 

0.8 


—^^—,—,—1—1—^— 

- 6-5 . 

-1 -,-1-^-,-1-1-^- 

‘ Serpens SMM 1 

: . J 

\ ^ 

_ kifVi'l.h I'M -.ibl-n ai1 nnJI ■ r^j ; 


snnrr ^pin-- 

^ Vela IRS 19‘ 


\ NGC 2071 (xO.4). 

L1641 S3 MMS 1 ' 


.. 

IRAS 20050+2720' 

rin„ n n. n 

LJV_J-U-“ - □ - - □ - UI-lI - ljulJ 

_ ^^^^^^^^^^^^^^^^_ 


-20 0 20 40 


Vlsr (km/s) 


^ RMS = Icr for channel width of 0.5 km s ^ 


sion at the off position. A spectrum taken at the off position re¬ 
vealed that the absorption is narrow and not affecting emission 
at the velocities of Vela IRS 19 which is offset by more than 10 
km s“^. 

3. Results 

3.1. Line Profile of central position 

Table presents the integrated intensities, peak temperatures 
and effective noise levels of spectra extracted from the central 


Fig. 1. CO 6-5 spectra taken at the central position. The baseline is 
shown in red. The line profile for LI641 S3 MMS 1 shows an example 
of the gaussian component fit (overplot in blue). The green lines show 
the source velocity, also listed in Table 1. 


(0,0) positions, assumed to be the gravitational centers. We note 
that not all protostars are covered by the beam at (0,0). E.g., 


NGC 2071-C ( jvan Kempen et al.|2012| ) and IRAS 20050-2720 


OVRO 2 (BOS) are located over 9" away from this position. Re¬ 
sulting spectra are shown in Figures (CO 6-5),|^(^^CO 6-5) 
and0([CI]). 


All lines are detected with S/N >10, with the exception of 
a non-detection of [Cl] in LI641 S3 MMS 1. Line profiles of 
the 6-5 are dominated by strong line wings, indicative of 
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Table 3. Parameters of the CO 6-5 component fits. 


Source 

^peak 

(K) 

Broad Component 
FWHM f TdV 
(kms“') (Kkms“^) 

^peak 

(K) 

Medium Component 
FWHM f TdV 
(kms“') (Kkms“/ 

^peak 

(K) 

Narrow Component 
FWHM f TdV 
(km s“^) (Kkms“') 

NGC 2071 

20.1 

27.9 

596.9 

51.3 

9.3 

508.5 

- 

- 

- 

L1641 S3 MMS 1 

2.8 

16.0 

47.5 

5.9 

5.6 

35.0 

- 

- 

- 

Vela IRS 17 

- 

- 

- 

28.3 

8.9 

267.8 

- 

- 

- 

Vela IRS 19 

- 

- 

- 

5.3 

13.1 

73.8 

2.5 

3.2 

8.5 

IRAS 20050-f2720 

- 

- 

- 

12.3 

12.1 

159.1 

- 

- 

- 

Serpens SMMl 

- 

- 

- 

9.3 

13.0 

129.1 

12.7 

4.6 

62.2 



Vlsr (km/s) 


Serpens SMM 1 




30 


[Cl] "P2-"Pi 


Vela IRS 17 




20 


10 - 





L164L S3 MMS 1. 



-20 


20 


40 


Vlsr (km/s) 


Fig. 2. 6-5 spectra taken at the central position. The baseline is 

shown in red. The green lines show the source velocity, also listed in 
Table 1. 


Fig. 3. [Cl] ^P 2 -^Pi spectra taken at the central position. The baseline 
is shown in red. The absorption feature at ~0 km s“^ in Vela IRS 19 is 
from the off-position but does not affect the main line emission. 


outflow activity. Absorption features are seen near the source 
velocities in all sources but their shapes differ from source to 
source. For NGC 2071, the absorption can be directly associated 
to large-scale material. The absorption is detected at three dis¬ 
tinct velocities coinciding with the velocities of the large-scale 
CO , and C^O 3 -2 emission profiles (12, 8 and 4 km s“\ 
see |Buckle et al.|20T0| ). As such, it is safe to assume absorptions 
are caused by cold material in the outer envelope and/or large- 
scale cloud. The and [Cl] lines are dominated by narrow 
emission components. Wider components in these lines are only 
seen for Vela IRS 19 (^^CO), NCC 2071 (^^CO and [Cl]) and 
Vela IRS 17 ([Cl]) 

To better analyze the different components, spectra are de¬ 
composed by fitting up to three Caussians to each profile: a 


‘narrow’ (FWHM <4 km s“^), a ‘medium’ (between 4 and 15 
km s“^) and a ‘broad’ (> 15 km s“^) component. Absorption 
features are corrected for. This method is similar to the meth¬ 
ods used by [Kristensen et al. p010b| ), [Kristensen et al.| p012[ ) 
and jSan Jose-Carcia et al. p()13| ) to deconstruct H 2 O and high-/ 
CO line profiles. It is possible to identify different outflow com¬ 
ponents using the ‘medium’ and ‘broad’ components and sepa¬ 
rate them from quiescent components traced by ‘narrow’ com¬ 
ponents. Results of the decomposition are presented in Table 
for CO 6-5 and Table 0 for ‘^CO 6-5 and [Cl] 2-1. A visual 
example is overplotted in Fig. [T]on the L1641 S3 MMS 1 spec- 
trum. 

For CO 6-5, ‘medium’ components dominate. Only 
NGC 2071 and L1641 S3 MMSl show ‘broad’ components, 
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Table 4. Parameters of the [Cl] 2-1 and 6-5 component fits. 


Source 

^peak 

(K) 

'^CO 
FWHM 
(km s“^) 

f 

(Kkms-b 

NGC 2071 

13.9 

8.0 

118.2 


9.1 

2.1 

20.2 

L1641 S3 MMS 1 

5.3 

1.8 

10.3 

Vela IRS 17 

15.1 

3.9 

62.8 

Vela IRS 19 

0.8 

8.7 

7.5 


3.2 

2.2 

7.4 

IRAS 20050+2720 

3.4 

2.9 

10.6 

Serpens SMMl 

0.9 

8.0 

8.0 


6.4 

2.6 

17.3 

Source 


[Cl] 2- 

1 


^peak 

FWHM 

jTdV 


(K) 

(km s“^) 

(K km s-‘) 

NGC 2071 

3.5 

10.9 

40.7 

L1641 S3 MMS 1 


no detection 

Vela IRS 17 

3.7 

7.6 

29.8 


5.0 

2.6 

13.5 

Vela IRS 19 

0.8 

6.5 

5.6 


3.4 

3.6 

13.3 

IRAS 20050+2720 

4.4 

2.7 

12.7 

Serpens SMM 1 

2.9 

3.2 

9.9 


while ‘narrow’ emission components were found for Vela IRS 19 
and Serpens SMMl. 

[Cl] and lines decompose into ‘narrow’ and ‘medium’ 
components, although the widths of the ‘medium’ components 
are often lower than corresponding ‘medium’ components. 
Observed variations in line width of ‘narrow’ components are 
caused by uncertainties in the fitting routine and the achieved 
S/N. 

3.2. Maps 

Fig. 1^ shows the CO 6-5 emission associated with the molec¬ 
ular outfiows in comparison the integrated emission within 2 
km s“^ of the Vlsr- Outfiow emission was measured by inte¬ 
grating between velocities of +4 to +20 km s“^ with respect 
to the source velocity. For the very broad fiow emerging from 
NGC 2071, these cuts were changed to ± 10 to ±40 km s“^ If 
known, positions of (sub)millimeter detected protostars are plot¬ 
ted with white crosses. For the Vela sources, no interferometric 
(sub)millimeter observations exist to identify individual proto¬ 
stars. Figure shows the emission of ^^CO (contours) and [Cl] 
(colors). Map sizes in these lines are typically smaller than the 
i^CO 6-5. 


Table 5. Derived ^^CO 6-5 average optical depths in the center and line 
wings. 


Source 

^center 

^blue 

^red 

NGC 2071 

>15 

2 

< 1 

L1641 S3 MMS 1 

>20 

2 

2 

Vela IRS 17 

>15 

4.5 

2.5 

Vela IRS 19 

>20 

- 

6.5 

IRAS 20050+2720 

>10.0 

<1.0 

<1.0 

Serpens SMMl 

>9.5 

2.2 

1.6 


compositions we conclude that the relative contribution of the 
‘broad’ component to the total integrated line fiux increases as a 
function of excitation energy. 

For sources in our sample where ‘broad’ components are de¬ 
tected in CO 10-9, counterparts in CO 6-5 are likely hidden by 
the noise; The focus for the CO 6-5 observations above was the 
size of the maps and not the sensitivity. 

^^CO emission is dominated by ‘narrow’ emission. As seen 
in the contours of Fig. it originates in circumstellar envelope 
in most cases. The exception appears to be IRAS 20050+2720, 
where ^^CO peaks 30" south-east of the protostars. For three 
sources, NGC 2071, Vela IRS 19 and Serpens SMM 1, a 
‘medium’ component is detected, but apart from NGC 2071, this 
component is much weaker than the ‘narrow’ component. 

The [Cl] emission is dominated by ‘narrow’ emission com¬ 
ponents. Only NGC 2071 shows emission with a FWHM >8 km 


4.2. Line Luminosities 


CO line luminosities show a relatively tight relation across the 
large range of luminosity and mass involved in star formation 
dWu et al.|[2005l [20T0l |San Jose-Garcia et al.|[MT^ . The line 
luminosity relation is defined as 


L = l(fL7 , 

bol 


( 1 ) 


with L the line luminosity (Lqo, or for CO 6-5, [Cl] 
^P 2 -^Pi and ^^CO 6-5 respectively. 

Figure shows the line luminosities compared t o Lboi, in- 
cluding results of CO 6-5 of low-mass protostars from |van Kem- 
|pen et al.] ( |2009b| ) and [van Kempen et al.| ( |2009a| ). An average 
va lue y of 0.83 +0.05 is f ound, almost identical to 0.84 +0.06 
of |San Jose-Garcia et al.| ( [20T3] ) using Herschel CO 10-9 line 
observations and spanning the full range in luminosity between 
1 and 10^ Lboi. We note that even with the small sample size, 
[Cl] line luminosities follow the expected correlation between 
low- and high-mass star formation through total luminosity with 
a slope of 0.96. 


4. Analysis 

4.1. Components at central position 

The decomposition of CO 6-5 in Tabl e [3] differs from both 
the CO 10-9 and 3-2 decompositions in |San Jose-Garcia et al 
(2013). ‘Broad’ components (> 20 km s ^ in width) dominate 
the CO 10-9, ‘medium’ components the CO 6-5, and ‘narrow’ 
components the CO 3-2. ‘Broad’ components are detected for 
CO 3-2, but are relatively weaker than the ‘narrow’ component. 
Similarly, the two detected ‘broad’ components are weaker than 
their CO 10-9 counterparts. From a comparison of all three de- 


4.3. Optical depths 

Table [^presents three different optical depths derived using the 
CO 6-5/^^CO 6-5 line ratios: at line center (reenter) and in each of 
the wings (Tbiue/^red)- A standard ^^CO/^^CO isotopologue ratio 
of 65 was used ( [Wilson & Rood|1994| ). 

Optical depths at line center range from 10 to >25, and are 
set by foreground absorption. In the wings, optical depths range 
between <1 and 6.5 with most values between 1 and 3. We note 
that these are upper limits as a result of the lack of signal at 
higher velocities in ^^CO 6-5. These optical depths are higher 
than upper limits for low-mass sources (^^1, e.g., NGC 1333: 
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Fig. 4. Maps showing the CO J-6-5 integrated emission in -20 to -4 {blue) and +A to +2Q {red) km s"^ bins to visualize the outflowing gas, 
overplotted on the integrated emission in a -3 to +3 km s"^ bin {dotted line and grayscale), representing the quiescent emission. All velocities are 
with respect to the individual source velocity (see Table 1). NGC 2071 was characterized by deflning outflow bins at -40 to -10 and -hlO to -h40 
km s"^ instead of the bins above. All contours, including those of the outflowing gas are normalized towards the peak intensity of the quiescent 
gas component at the central position (Tpeak in Table 2). Levels are in turn given in 10%, 20%, ..., 80%, 90% w.r.t. to this peak intensity. Where 
known, locations of (sub)millimeter interferometry sources are shown with ‘x’. 
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Fig. 5. Integrated intensity of 6-5 (contours) overplotted on the integrated intensity of [Cl] (colorscale). Both distributions are normalized 
to the peak integrated intensity of that particular line. For the contours are in levels of 10%, 20%,..., 80%, 90% w.r.t. to this peak intensity, 
with the lowest contour higher than 3 times the noise level in Table 2. The [Cl] emission scales between 3 times the noise level and the highest 
intensity in the map, which can be found in Table 2. Where known, locations of (sub)millimeter interferometry sources are shown with ‘x’, except 
for the Vela sources. 
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Table 6. Line wing ratios of 6-5!J 2 - J\ and temperature and density estimates using RADEX. The density estimate is the lower limit for 

thermalized emission. The temperature assumes a density n of 10^ cm“^. 


Source Name 

J 2 -J 1 

Ref.‘ 

Ratio 

T 

(K) 

Blue 
n(H2) 
(10^ cm"^) 

Ratio 

T 

(K) 

Red 

n(H2) 

(10^ cm“^) 

NGC 2071 

3-2 

1,2 

1.0 

> 100 

300 

1.2 

>100 

100 

L1641 S3 MMS 1 

3-2 

2 

2 

>70 

500 

2.0 

>70 

400 

Vela IRS 17 

1-0 

3 

1.5 

90 

9 

1 

>100 

9 


1-0 

4 

<1 

100 

10 

<0.8 

>140 

25 

Vela IRS 19 

1-0 

3 

<1 

>80 

10 

1.4 

75 

10 


1-0 

4 

<1 

>80 

10 

1.4 

75 

10 

IRAS 20050+2720 

2-1^ 

5 

1.5 

>70 

400 

3 

>50 

200 


2-1^3 

5 

>10 

>10 

10 

>4 

>50 

75 

Serpens SMMl 

4-3 

6 

1.7 

>80 

500 

3 

>60 

100 


3-2 

2,7,8 

1.0 

>90 

300 

1.0 

>50 

<300 


3-2^ 

2,7,8 

2.5 

>60 

100 

- 

- 

- 


2-1 

9,10 

0.9 

>80 

200 

1.5 

>70 

<200 


^ l:|Buckle et al.|j2010|, 2:|San Jose-Garcia et al. (|2013|, 3:|Wouterloot & Brand|Jl999|, 4:|Elia et al.|j2007|, 5:|Bachiller et al.Hl995|, 6:|Hogerheijde et al.|(|l999|, 
7: |Dionatos et~ni(l2^. s'lGraves et al.ld2010> , 9 lWhite et al.K1995h 10: [Davis etarifT99^ 

^ at off-positions of |Bachiller et al.|jl995) . 


2 

^ EHV component at identified positions 
^ Tkin? assuming nH 2 <10^ cm“^. 


Yildiz et al.|2Q12 HH46: van Kempen et al.|2009b| ). Exceptions 


are IRAS 20050+2720 and NGC 2071, which show optically 
thin emission in the line wings. 

Optical depths at different positions are consistently equal or 
lower than the optical depths given in Table For convenience, 
the optical depth for the flows are used for the full maps. 


4.4. Outflow properties 

Outflow parameters such as kinetic temperatures, densities, out¬ 
flow forces and kinetic luminosities can be derived by calcu¬ 
lating the non-LTE radiative transfer parameters and comparing 
those with observed line emission. For this sample, the off-line 
version of the RADEX code ( |van der Tak et al.|2007| ) was used 
provide constraints using ratios of the observed CO 6-5 emission 
over previously observed transitions of lower excitation (see Ta¬ 
ble [^. Diagnostical plots for ratios with respect to CO 6-5 are 
presented in Fig. [ 7 ] 


4.4.1. Temperature and density 


Diagnostical plots produced by RADEX as shown in |van der Tak| 


|et al.| ( |2007| ) reveal that ratios provide solutions with a degener¬ 
acy between temperature and density. To break this degeneracy 
other molecular tracers are required to independently derive ex¬ 
citation constraints. CO line ratios covering three or more tran¬ 
sitions can provide additional information, but are often insuffi¬ 
cient to completely solve the degeneracy. However, ratios using 
the CO 6-5 line emission are able to exclude large areas of the 
parameter space. E.g., tempe ratures under 50 K are found to be 
excluded for many outflows ( |van Kempen et al.|2009b|at |Yildiz 


et al.||2012| ). For a more thorough discussion on RADEX solu¬ 


tions concerning CO and including thermal versus sub-thermal 
excitation, we refer the reader to |Yildiz et al.| ( |2012| ). 

To derive the excitation parameters of this sample, RADEX 
was run in the optically thin limit, adopting the following param¬ 
eters: a fine width of 10 km s“\ a column density of 10^^ cm“^, 
and a background radiation field of 2.73 K. 
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In turn, the excitation conditions were investigated by con¬ 
sidering two scenarios. First, (lower limits to) the temperatures 
were derived by assuming a density of 10^ cm“^. Second, a lower 
limit for the density is given. This solution is the lowest density 
at which emission is fully thermalized. In other words, the den¬ 
sity given is the lowest density for which the ratio solely depends 
on temperature (the limits found at the right sides of the diagnos¬ 
tical plots in Fig. [ 7 ] where lines are horizontal). These two sce¬ 
narios were chosen as they likely best reflect the true physical 
conditions. Results can be found in Table (6] 

Temperature limits of 50 K are found for all flow^ These are 
in ag reement with temperature constraints for low -mass proto¬ 
stars ( [van Kempen et al.|2009b|at|Yildiz et al.|2012| ) Technically, 
lower temperatures are not excluded, but require densities > 10^ 
cm“^. Spherical envelope modelling restricts such densit ies to 
<2000 AU from the central protostar (See Dusty models of |Kris- 
tensen et ar]|2012| ). Similarly, BOS restricted densities in IRAS 
20050+2720 to 1.3 - 3x 10^ cm“^ at radii <3000 AU. In general, 
densities at larger radii are significantly lower, with most cloud 
densities derived to be on the order of lO'^ cm“^. Compression 
factors can be invoked to compensate for this change in density 
(i.e, local density enhancements). However, compression factors 
of more than three orders of magnitude are required to keep CO 
emission thermalized along the entire observed flow (>20,000 
AU). Molecular tracers such as H 2 O ( [Santangelo et al.|2012| ) or 
HCO^ (van Kempen et al.|2009b| ) independently provide density 
constraints of <10^ cm“'^, limiting compression factors to two 
orders of magnitude. As such, the sub-thermal low-^, high-Y so¬ 
lution is the preferred solution. 

4.4.2. The optically thin limit 

Tests were carried out to verify the optically thin assumption 
used above and the effect higher optical depths would have on 


^ The limit of 10 K for the blue flow of IRAS 20050 was derived with 
a synthesized beam and is likely suffering from filtered out large-scale 
emission. 
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^bol [L©] 



^bol [L©] 






Fig. 7. RADEX diagnostic plots of CO 2-1 {top), 3-2 {middle) and 4- 
3 {bottom). The derived line ratios of the blue flow of Serpens SMMl 
are highlighted in red to serve as an illustration for the sub-thermal and 
thermal excitation scenarios. It is evident that temperatures below 50 K 
are excluded in both limits as ratios of CO 2-1 and 3-2 over 6-5 do not 
agree with those of 4-3. 


Fig. 6. Top: CO 6-5 line luminosity vs Lboi, Middle: ^^CO 6-5 line 
luminosity vs Lboi, Bottom: [Cl] ^P 2 -^Pi line luminosity vs Lboi- Fits to 
the line luminosities are shown in red, with the slope labelled ’y’. The 
slopes are within the error bars of |San Jose-Garcia et al.|(|2013| ). Open 
symbols are the low-m ass sample fr om Yildiz et al.| ( |2013| ). Note that 
no [Cl] was reported i n |Yildiz et al.| ( [T013| ). The [Cl] detections of |van| 
|Kempen et al.| ( |200^ are included. 

the density and temperature derivations. This was done by us¬ 
ing significantly higher column densities cm“^) in 


the RADEX simulations. These optically thick solutions provide 
significantly higher constraints on temperature (>200 K instead 
of >50 K), while the effective ^crit was found to increase, thus 
increasing the lower limit on the density given in Table Since 
the derived optical depths for outflowing gas are mostly upper 
limits, adopting the optically thin RADEX solutions provides us 
with the most conservative, but likely more realistic, estimate. 
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Tabl e 8. Correction factors use d to multiply the uncorrected value with, 


from 

text. 


Downes & Cabrit 


( |2007| ^ for an inclination i of 30 degrees. See 


u 

Fco 

3cr 

0.29 

2.S 

1.4 


N= 10' 


3k^ 


T + 


hv 


^ 6/f(// + l) 






^-HvIkT 


f 


Tmb 


-dV 


Fco = cx 


K(Y,i[fT^bTi=^V'dV’]iWm 




lobe 


1 . 


Densities derived from line ratios were found to be 10"^ 
or higher. Therefore, the 77 = 1 case of 


cm 


Downes & 


Cabrit] ( [2007 1 ) is used instead of the geometric mean of 0.1 


and 


1 adopted by [van der Marel et J ( |2013| ). The latter 


i = 90-a, as used in |Downes & Cabrit] ( |2007) . 


4.4.3. Velocity and spatial variations 

Line ratios were found to be relatively constant in velocity, with 
ratio variations typically on the order of 20 to 40%. This is simi¬ 
lar to the behaviour of CO in flows around HH 46 and NGC 1333 
IRAS 2 (Fig. 10 in |van Kempen et al.|2009b[|Yildiz et al.|2012| ). 

Line ratios at other positions also show little to no signifi¬ 
cant difference. It can thus be concluded that variations in the 
excitation mechanisms along the large-scale flows are small. 


4.4.4. Mass outflow rate and outflow force 

Outflow forces and kinetic luminosities are derived using the 
H 2 column density. Column densities are derived using Eq. 1 
of |Hogerheijde et al.| ( |1998| ) using parameters for the CO 6-5 
transition: 


( 2 ) 


where k is the Boltzmann constant, h the Planck constant, p 
the permanent dipole moment (0.122 Debye for CO), v the fre¬ 
quency of the transition, Ji the quantum number of the lower 
rotational state, and f - e~'^))dV the integrated line in¬ 

tensity, corrected for optical depth. All quantities are in cgs units, 
except for the velocity which is in km s“^ The total mass is cal¬ 
culated by summing column densities across the map, assuming 
a H 2 /CO ratio of lO^. 

The outflow force, Fco^ is derived using the integrated in¬ 
tensity as a function of velocity, corrected for optical depth and 
subsequently integrated over the observed area of pixels i and 
in turn corrected for the inclination. From recent benchmark¬ 
ing, the most reliable method to calculate outflow forces is the 
’M7’ or ’separation’ method ( [van der Marel et al.|2Q13| ). In this 
method the dynamical age and force of a flow are considered 
to be independent quantities. The dynamical age, is defined 
as the measured radius, R, divided by Vmax- Using the intensity 
weighted velocities the outflow force is thus expressed with the 
following equation: 


corresponds to densities below 10"^ cm“'^. As before, these 
low densities would imply unlikely temperatures of 200 K 
or higher. Although no t excluded, typical envelope models 
( [Kristensen et al.|2012| ) already have higher densities. In ad¬ 
dition, extrapolation of the correction factors from [Downes 
|& Cabrit] ( 120071 ) above 100 K is not reliable. 

2 . AVniax was measured using a 3 (T limit, instead of a Icr limit. 
Data quality at 690 GHz was found to be insufficient to make 
a reliable Icr limit derivation. With the higher system temper¬ 
atures of the CHAMP^ due to the lower atmospheric trans¬ 
mission at 690 GHz, the effective S/N of the CO 6-5 obser¬ 
vations here are a factor of 5 or more lower than for CO 3-2 
used by ] van der Marel et al.] ( ]2013] ). To avoid any potential 
systematic errors introduced by the data quality but still cor¬ 
rectly approach true values for A Vmax correctly, both and 
Fco were corrected with an additional factor of 1.4. This was 
tested by extrapolating the gaussian fits to Vmax and found to 
be robust. The correction factor was derived from tests us¬ 
ing the appendix of ] van der Marel et al.](]2013]) as well as 
similar tests on data presented here and van Kempen et al. 
P009b] ). td is divided by this correction factor, while Fqo is 
multiplied. 

3. No reliable information on individual viewing angles of the 
outflows with respect to the plane of the sky is available. An 
average value of 32 degrees for the angle of the outflow with 
the plane of the sky is adopted. This is the expected mean 
value for a randomly distributed sample of outflow inclina¬ 
tions. 


Table lists the final correction factors used. As a comple¬ 
ment to the outflow force, the kinetic luminosity of the flows, 
^kin, was calculated using 


^kin — ^CO ^ Fmax/2. 


(4) 


Table |7] lists the final values for all outflow parameters. Vmax 
values of both lobes are consistently of similar value, with the 
exception of the blue side of Vela IRS 19, which is not detected 
strongly. Its derived parameters are considered lower limits. Dy¬ 
namical times are a factor of 2 shorter than the average for the 
low-mass protosta r sample (van Kempen et al.]]2009ct ]van der 
Marel et al.|2Q13] ). This may be indicative that no evolved inter¬ 
mediate mass protostellar cluster was included. Outflow forces 
are factors of 10 to 300 higher than low-mass protostars (Bon- 


temps et^ 


19961. When compared to the results of Duarte 


( 3 ) 


Cabral et aL](2013] ), the outflow forces found are consistent with 
the low-end of values for high-mass sources (10“'^ M© yr“^ 


km 


Here c is the inclination correction and K the temperature- 
dependent correction factor. Riobe is the radius of the lobe. For 
more information on this method, see ]van der Marel et al.](201^ . 
The correction factor is derived from the values of Table 6 of 
[Downes & Cabrit] (2007] ) (see Table 7). 

Owing to atmospheric effects, observations at 691 GHz can¬ 
not be obtained with a similar signal to noise ratio as its low-/ 
counterparts within reasonable times. As such, the M7 method 
was changed on three points w.r.t.]van der Marel et al.]([20T^: 


s ^ or higher for luminosities of 100 and higher). 

4.4.5. The multiple flows of IRAS 20050+2720 

BOS observed IRAS 20050+2720 using a resolution of 3 arcsec- 
onds. This allows one to identify all individual outflows. It is 
seen that the observations of interferometers do not resolve out 
any emission. The total outflowing mass and outflow force co¬ 
added for both flows in IRAS 20050+2720 add up to 0.25 M© 
and 6.4 10“"^ M© km s“^ yr“^ (BOS). Co-added outflow forces 
from CO 6-5 are less than a factor of 3 higher. With the sys¬ 
tematic difference due to the difference in temperature (29 K for 
BOS, 50 K for this study), this factor is well within the assumed 
inaccuracy of the M7 method (van der Marel et al.|2013] ). 
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Fig. 8. Luminosity versus outflow force of the blue (circle) and red (square) lobes. Black symbols are the derived values for the intermediate 
mass protosta rs, dark grey symbols are v alues from |Yildiz et al.|f2015| ), li ght grey symbols are value s from |van Kempen et al.| ( |200^ and white 
symbols from |van der Marel et al.| ( [20T3] ). All values except the ones from van der Marel et aL](|2013| ) are derived using CO 6-5. White diamonds 
are the values from BOS. The red line represents the relation proposed by |Bontemps et al.| ( |1996| ). 



^env [M©] 


Fig. 9. Envelope masses versus outflow force of the blue (circle) and red (square) lobes. Black symbols are t he derived values for the in termediate 
mass protostars, dark grey symbols are values from |Yildiz et al.| ( f2015| ), light grey symbols are values from |van Kempen et al.| ( |200^ and white 
symbols from|van der Marel et al.|(|2013|). White diamonds are the values from BOS. The red line represents the relation proposed by|Bontemps| 
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Table 7. Outflow Parameters. The dynamical time, td, and outflow force, Fco were corrected for inclination using the factors listed in Table 9 as 
well as a factor of 1.4 to compensate the method in determining Vmax- See text. 


Source 

Mass 

Umax 

R 

td 

M 

Fco 

^kin 


(Mo) 

(km s“^) 

(10^ AU) (10^ yr) 

(10-5 Mo 

(10-4 Mo 

(Lo) 






yr-') 

yr"^ km s"^) 


Red Lobes 

NGC 2071 

0.39 

39.0 

4.2 

1.1 

36.3 

145 

45.8 

L1641 S3 MMS 1 

0.022 

12.0 

0.93 

0.77 

2.8 

3.5 

0.34 

Vela IRS 17 

0.27 

10.0 

3.5 

3.5 

7.7 

7.9 

0.64 

Vela IRS 19 

0.38 

13.0 

2.1 

1.6 

23.8 

31.9 

3.4 

IRAS 20050-^2720^ 

0.31 

14.0 

4.9 

3.5 

8.9 

12.8 

1.5 

Serpens SMMl 

0.006 

12.5 

1.3 

1.1 

0.5 

0.6 

0.06 

Blue Lobes 

NGC 2071 

0.41 

36.5 

42 

1.2 

35.3 

130 

38.4 

L1641 S3 MMS 1 

0.017 

14.5 

27 

1.9 

0.9 

1.4 

0.16 

Vela IRS 17 

0.31 

10.0 

42 

4.2 

4.6 

4.6 

0.37 

Vela IRS 19^ 

0.05 

>3.5 

27 

<7.7 

0.7 

>0.23 

>0.007 

IRAS 20050-f2720i 

0.22 

14.0 

42 

3.0 

7.4 

10.6 

1.2 

Serpens SMMl 

0.005 

16.5 

9 

0.5 

1.5 

2.6 

0.35 


^ average over the detected flows in |Bachiller et al.|fl995) and BOS. 
^ strongly limited by detection of Vmax- 


5. Discussion 


5.1. Does fragmentation enhance outflow forces? 

Using a large sample o f outflows emerging fro m low-mass pro- 
tostellar environments, [Bontemps et J^(|1996| from here on re¬ 


ferred to as B96) revealed a relation between the bolometric lu¬ 
minosity of the driving source and the force of its outflow (from 
here on referred to as the B96 relation). The relation inherently 
has a relatively large scatter, and the influence of evolutionary ef¬ 
fects could not be determined accurately. Since then, many other 
studies corroborated this relation, and reveals it may apply across 
the many orders of magnitude in luminosity in star formation, up 


to and including high-mass star formation (Duarte-Cabral et al. 
2Q13| ). Figure shows the B96 relation (red line) in compari¬ 


son with the results obtained here (black fllled points) and BOS 
(white diamonds). Studies of low-mass protostars using CO 6- 
5 and/or similar methodology to derive outflow forces are in- 
cluded as reference (|van Kempen et ar]|2009a[ |van der Marel 


et al.||2Q13HYildiz et al.||2015[ grey S 3 mbols). For the interme¬ 


diate mass flows, the outflow force at first sight correlates with 
the total Lboi with the B96 relation, similar to the conclusions of 
BOS. The largest deviation from the relation is the NGC 2071 
flow, where the outflow force is an order of magnitude higher 
than expected based on the B96 relation but still well within the 
observed scatter. 

Interestingly enough, the median of the the outflow forces 
in [Yildiz et al.| ( |2015| ) as derived for low-mass sources using 
CO 6-5 is also statistically significant above the B96 relation. 
An excess of an order of magnitude for fragmented intermediate 
mass sources was identified earlier by BOS, although their results 
were inferred from data with a higher uncertainty on the outflow 
forces. Direct calculations of gas temperatures, correct deriva¬ 
tions of optical depths and a uniform derivation of Fco have im¬ 
proved the reliability of the derived values. The observed excess 
of BOS can thus neither be corroborated or invalidated. What is 
more likely is that with the improvements made in observations, 
the better understanding of the derivation of outflow forces and 
the access to tracer lines better suited to track entrained outflow 
material (in this case the CO 6-5), outflow forces are slightly 
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higher than the original B96 relation, although the most impor¬ 
tant aspect of it, the slope, remains the same. 

Fi gure [9| compares the Fqo and Menv relation of |Bontem^ 
et al.| ( |1996| ) for the same set of samples above with the B96 
relation shown as a red line. The correlation between Mgnv and 
Fco is clearly visible for intermediate mass sources and in direct 
agreement with those derived from low-mass sources. It should 
be noted that compared to the low-mass sources (shown in white, 
light gray and dark gray), the scatter for the intermediate mass 
sources (shown in black or white diamonds) has increased by 
a factor 2. Most likely this enhanced scatter is caused by mea¬ 
surement uncertainties due to intermediate mass sources being 
more distant. Intermediate mass sources are a factor of 2 to 5 
more distant than typical low-mass sources. From these results, 
the results of BOS cannot be corroborated. The observed en¬ 
hancement in outflow forces seen by BOS are reproduced, but 
are within the scatter of the B96 relation. In addition, the ad¬ 
vancement of more accurate observations and differences be¬ 
tween them |van der Marel et ak] ( |2013| )and corresponding con¬ 
straints can be invoked to explain any changes . 


5.2. Neutral carbon 


Neutral carbon has long been assumed to be created from inter¬ 
actions of the gas with the ISRF (Interstellar Radiation Field), 
which produces atomic gas components in a PDR scenario 
from photo-dissociation of carbon-bearing species ( [Hollenba^ 
& Tielens|1997| ), although atomic gas within the outflow c avity 


surface PDRs could contribute to neutral carbon emission ( |Hol-| 
[lenbach & Tielens|1997| ). 

Indeed, the observed spatial distribution of the [Cl] ^P 2 - 
^Pi emission is clearly different from that of CO. In addition, 
no correlation to the direction and/or strength of the outflow is 
detected. Figure revealed [Cl] emission to be smoothly dis¬ 
tributed over the circumcluster envelope, with concentrations at 
or near the protostellar positions. The only correlation seen for 
[Cl] is in the line luminosity, which reproduces the same slope 
as derived for CO and 6-5. This is not surprising as [Cl] 
is clearly expected to be coupled to CO, which scales linearly 
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with luminosity. Whether or not this correlation is evidence for 
atomic gas emission being correlated with outflow activity itself 
cannot be confirmed owing to the low number of detections of 
[Cl]. 


6. Conclusions 

This paper presents new spectral line observations of six pro- 
tostellar clusters of intermediate mass, with total luminosities 
ranging from 30 L© to ^^750 Lq. CO /=6-5, J=6-5 and 

[Cl] spectrally resolved maps were obtained with the CHAMP^ 
instrument on APEX. Using line decomposition, accurate opti¬ 
cal depths and line luminosity relations, densities, temperatures, 
forces and kinetic luminosities of the molecular outflows were 
derived and presented. The conclusions can be summed up as 
follows: 


The CO 6-5 line profiles are dominated by outflow related 
emission, but show quiescent emission in the cloud as well. 
Mid-/ CO line luminosities adhere to the correlation be- 


tween total luminosity a nd line luminosity identified by jSan 
Jose-Garcia et al.| ( [20131 ) for low- and high-/ CO. 


There is no corroboration of the result presented in Beltran 
|et al.| ( |2008| ) that proposed an apparent enhancement in out 
flow force for fragmented intermediate mass sources. Al¬ 
though an enhancement of the outflow forces as a function 
of total bolometric luminosity is seen in comparison with the 
original B96 relation, this increase can also be attributed to 
methodology or the improvement in temperature and density 
derivations due to the inclusion of mid-/ CO. 


Future work on outflows emerging from protostellar clusters of 
intermediate mass require observations down to scales of indi¬ 
vidual protostars. Properties of individual sources are necessary 
to draw conclusions about the influence of fragmentation. Al¬ 
though near and mid-infrared observations have been acquired 
that can be used (e.g., Spitzer, WISE), (sub)millimeter obser¬ 
vations with sufficient spatial and spectral resolution are rare. 
ALMA is able to routinely do such observations in minutes 
through several CO transitions, although other interferometers at 
these wavelengths (SMA, CARMA, IRAM Plateau de Bure and 
its successor NOEMA) should not be discounted even though 
their lack of access to mid-/ CO lines will limit their effective¬ 
ness. The GREAT instrument on SOFIA and in particular the up- 
GREAT array extension of the instrument may spectrally resolve 
foreground atomic gas from outflowing atomic gas. In combina¬ 
tion with ALMA Band 8 and 10, such observations must be used 
to interpret the [Cl] observations. 
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